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We present an exact real-space renormalization group (RSRG) method for evaluating the dy- 
namic structure factor of an infinite one- dimensional quasiperiodic period-doubling (PD) lattice. 
We observe that for every normal mode frequency of the chain, the dynamic structure factor S(q, uf) 
always exhibits periodicity with respect to the wave vector q and the presence of such periodicity 
even in absence of translational invariance in the system is quite surprising. Our analysis shows that 
this periodicity in S(q, u>) actually indicates the presence of delocalized phonon modes in the PD 
chain. The Brillouin Zones of the lattice are found to have a hierarchical structure and the dispersion 
relation gives both the acoustic as well as optical branches. The phonon dispersion curves have a 
nested structure and we have shown that it is actually the superposition of the dispersion curves of 
an infinite set of periodic lattices. 

PACS Numbers: 71.23.Ft, 61.44.-n, 63.20.Dj 



I. INTRODUCTION 

The discovery of icosahedral phase in Al-Mn alloys by Shechtman et al. 0] initiated experimental and theoretical 
research on quasiperiodic systems. Theoretical investigations are mostly based on one-dimensional quasiperiodic 
systems, as they can be described by simple models for pointing out the relevant properties of these systems. There 
are many interesting works in the literature concerning the electronic [|]^| arL d vibrational properties JlO| — p~5[l °f 
these lattices. At present there is considerable amount of understanding about the nature of the eigenfunctions 
and the eigenvalue spectrum for this class of deterministic hierarchical systems lacking translational invariance. The 
general consensus is that the properties like critical states, Cantor-set spectrum, power-law behavior of the density of 
states etc. are the characteristic features of the quasiperiodic or aperiodic systems, though there are evidences 
of extended states |B0 in such systems. Moreover, there are also some reports on the anomalous transport Q and 
optical localization ^] properties of these systems. But so far very less attention has been given regarding the dynamic 
response function of these lattices, a quantity which is directly connected to the scattering cross-section obtained from 
neutron diffraction or X-ray scattering experiments. Even the computation of the dynamic structure factor for one- 
dimensional quasiperiodic system is quite meaningful, as it is now possible to fabricate semiconducting superlattices 
epitaxially p.6fl , and one can perform neutron diffraction or X-ray scattering experiments on these samples. Therefore, 
we can directly compare the theoretical predictions with the experimental results. As we will see, the study of the 
dynamic structure factor also gives information about the nature of the eigenmodes and consequently one can predict 
the character of the modes from experimental observations. 

In this paper we have addressed the dynamical properties of quasiperiodic systems. We have computed the dy- 
namical structure factor for the quasiperiodic period-doubling (PD) chain using a technique based on the real-space 
renormalization group (RSRG) method. Here we present a formalism for calculating the dynamic structure factor of 
the system using ideas developed in Ref. fl^| . This is an exact method for computing the dynamic structure factor 
for a class of one-dimensional quasiperiodic chains, where the original lattice can be split into a finite number of 
self-similar sublattices. 

Our study of the dynamic structure factor also gives information about the nature of the eigenmodes. The dispersion 
relation for the PD chain can be easily obtained from the dynamic structure factor, and the normal mode frequencies 
become readily available. Hence it is possible to characterize the nature of the states using the transfer matrix 
technique. The most interesting observation is that for every normal mode frequency of the PD chain, the dynamic 
structure factor becomes a periodic function of the wavevector, and it actually signifies that this lattice supports 
delocalized classical vibrational modes for all these eigenmodes. We have also shown that for a certain choice of the 
parameters of the system, the dynamics of the PD chain exactly coincide with that of various periodic chains having 
different periods, and the corresponding features are reflected in the dispersion relation. Our study also reveals that 
one can infer about the delocalized nature of the states from a mere inspection of the X-ray or neutron scattering 
experiment data. 

In this work, though we have concentrated our attention on the quasiperiodic PD chain, the ideas are quite general 
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and will be applicable in many other quasiperiodic systems. We have studied the PD chain as a prototype example for 
illustrating our method. An added interest in taking the PD lattice is that the problem of the existence of delocalized 
modes in this lattice can also be studied on the basis of dimmer-type correlation |5|,[l7j among the atoms. But we 
would like to stress that though the method of Ref. || is theoretically very elegant, from the practical point of view its 
applicability is highly limited. The procedure of Ref. J5| actually gives only few delocalized normal mode frequencies, 
and the determination of the entire spectrum becomes a formidable computational task as it requires the solutions 
of polynomial equations with progressively higher and higher degree. On the other hand, the dynamic structure 
factor automatically gives the entire spectrum and the nature of the eigenmode corresponding to every normal mode 
frequency can be easily characterized. An interesting feature of the PD chain is that the trace map associated with 
this lattice leads to a polynomial invariant which should normally ensure the critical nature of all the eigenstates as 
well as a Cantor-set spectrum |3|,[fo]]. However, in this work we find that this lattice supports an infinity of delocalized 
states in spite of the Cantor-set nature of the spectrum. 

This paper has been organized as follows. We first describe our model in Section II. We present in Section III the 
RSRG scheme for finding the dynamic structure factor of a quasiperiodic system, while we briefly outline the transfer 
matrix method in the next Section. Section V is devoted for the analysis of our results and we conclude in Section 
VI. 



II. DESCRIPTION OF THE MODEL 

We describe the vibrational properties of the PD chain considering nearest neighbor harmonic interaction among 
the atoms. A portion of the PD chain is shown in Fig. 1. The PD sequence can be generated from two symbols L 
and S using recursively the substitution rule L — > LS and S — > LL. Thus L, LS, LSLL, LSLLLSLS, . . . etc. are 
the first few generations of the PD sequence. Now we can construct a PD chain by considering the symbols L and 
S as representing 'long' (L) and 'short' (S) bonds in the chain. In our model m a , mp and m 7 are the masses of the 
atoms flanked by L — L, L — S and S — L bonds respectively, and, fci and ks are respectively the spring constants 
across the long and short bonds. It should be noted that this kind of labeling of the sites and the bonds is essential 
for implementing the RSRG decimation procedure of Southern et al. |fl8|| . We shall see that successive renormalized 
lattices can be represented in terms of only these five parameters and the number of the parameters do not grow with 
iteration. In other words, within this labeling scheme, the form of the equations of motion of the system remains 
invariant under renormalization. The well-known on-site model and the bond-model for the PD chain turn out to be 
the special cases of this general situation. 

The expression for the dynamic structure factor S(q, lo) corresponding to the wavevector q and frequency u> can be 
written in terms of the single particle Green's functions as follows 

S(q,to) — lim lim lm(GN(q, oj + iS)) , (1) 

where Gjv(<Z,w) = (l/-^0 Eij & lq ^ ri ~ rj 'Gij . Here ri is the position vector of the ith atom and N represents the total 
number of atoms in the chain. In the harmonic approximation with nearest neighbor coupling between the atoms, 
the single particle Green's functions Gjj's satisfy the following equations of motion 

tiGij = —5ij + fci^+iGi+y + ki^-iGi-ij t (2) 

where = k i i+ i + fcj^—i — m^w 2 , m, being the mass of the ith atom, and fc™ the spring constant between the ith 
and jth atoms. 

Now the calculation of S(q, u>) for a translationally invariant system is very straight forward. Gjv(<7, u>) can be easily 
computed from a Fourier transformed version of the Eqs. (j^). However, for quasiperiodic systems the determination 
of the sum in Gjsr(q,u>) is not trivial. We shall show in the next section that this sum can be evaluated in a very 
elegant way using RSRG method. The calculation of S(q,u>) essentially amounts to the iteration of certain recursion 
relations and it can be computed with arbitrarily high accuracy. 

III. RSRG SCHEME FOR THE DETERMINATION OF S(Q,u) 

Let us now consider the sum G l {q,u) = ^.e^-^Gy. This sum is independent of the index i in a periodic 
system due to the translational invariance. However, in the case of a quasiperiodic system, it depends on the index i 
as no two sites are equivalent in the quasiperiodic lattice and from Eqs. (0) we obtain 
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e^fau;) = -F l + fe M+1 e-^ 1 - r -)G l+ H9,^) + A :M _ ie ^- r -- 1 )G^ 1 ( g , W ) , (3) 

where all F^s are initially equal to unity. The use of new notations Fi's for representing unity does not directly follow 
from Eqs. (|^). We have introduced these notations by hand and we will see how it facilitates the determination of 
S(q,co) in quasiperiodic lattices. Even though all i^'s are initially the same, they become different upon renormaliza- 
tion. However, from the symmetry of the lattice, we observe that there will be only three distinct type of fj's, and 
we can identify them as F a , Fp and F 7 corresponding to the a, (3 and 7 site of the lattice. 

It is now necessary to determine all G l (q,o;)'s from Eqs. (||) for evaluating Gjv(q, ui) since we have 

G N {q,u)^{l/N)Y,GKq^) . (4) 

i 

In the quasiperiodic PD chain G l (q,ujys are all distinct and instead of determining them directly from Eqs. (||), we 
use RSRG technique for finding them. We split the original chain into two self-similar sublattices SI and T as shown 
in Fig. 1. The fi-sublattice is obtained by eliminating the sites using the decimation rules LS — ► L and LL — » S, 
while the corresponding rules for the T-sublattice are SL — > L and LL — > S. Both f2 and T sublattices again form 
two new PD chains at some inflated length scale. All the sites of the original PD lattice are distributed among these 
two sublattices, and thus they are complementary to each other (see Fig. 1). This complementary nature of f2 and 
r sublattices also ensures that no information is lost by introducing this splitting procedure. It is now possible to 
generate two sets of renormalized equations for G z (g,w)'s, one corresponding to f2-sublattice while the other for the 
T-sublattice. The equations for fi-sublattice can be obtained from Eqs. (||) by eliminating all G*(c7,w)'s belonging to 
the T-sublattice, and the resulting equations can be cast in the same form as that of the original set of equations (^), 
provided we rename the sites appropriately and renormalize the parameters as follows 

e' a = e^-u (kl+k 2 s ), 

6 fs = € 1 - (klu) a + kgUp), 

e 7 = e a - k\{uj a +wp), 
F' a =F 7 + Fpk s cope lqas + FfikLLOpe'^ , 
F' = F^ + F a k L cu a e~ lq ^ + F p k s wpe iqas , 
F^=F a + F a k L Lo a e lqaL + FpkLiope- 1 ^ , 
k L = k L k S 0Jp , 

kg = kj{O a , 

a L = a L + a s , 

a s = a L + a,L ■ (5) 
Similarly for the T-sublattice the recursion relations can be written as 
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Here e, = fci,i+i + fci,i-i — rn,iU> 2 , u>i — 1/ej and i refers to a, (3 or 7. In Eqs. (J5J) and (^) and as respectively 
represent the 'long' and 'short' bond lengths. Let us represent the above two transformations by Tfi(p) an d denote 
the renormalized Green's functions as GL r j. Now we can recast GAr(q,w) as 
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G N {q,u) = p (n) G ( "\q,uj) +p { r)G ( P(q,u), (7) 

where G^\q,to) = (1/N\) J2iex G^iq.uj). Here A can be either Q or T, and is the number of sites in the 

renormalized tt(T) sublattice. The coefficients P(q\ and denote the fractions of total sites which belong to 

and r sublattices respectively, i.e., p(p\ = Nq/N and = Nr/N. The expressions for Gat((?,o;) and G^\q,uj) 
are structurally same, the former being defined in terms of the parameters of the original PD lattice while the latter 
refers to those of the renormalized A-sublattice. Eq. (|?|) shows that GAr(g,o;) is a linear combination of G^'(q, u) and 

G N (q, to) with coefficients p^ and P(r)j where p^ = p^) = 1/2- Since each of these sublattices again forms a new 
PD chain, we can treat them at the same footing as the original chain. Thus in the renormalized 0,(T) sublattice, 

G^ 1 ^ takes the same role as that of G^iq, oj) for the original lattice. As we can split further each of these new chains 

into f2 and T sub-sublattices, it is again possible to express both G^^w) and G^\q,uj) in the form of Eq. (Q). 
Thus we have 

G N {q,uj) = P(nn)G ( " n \q,uj) +P(nr)G^ r) (q, w) + p {rn) G^ n \q,u) + p ( rr)G ( £ r \q,Lj). (8) 

Here we denote the two branches resulting from fl sublattice as (00, OT) , while those from T-sublattice as (rf2,rr) 
. The coefficients can be written as pi^ v ) = with /i, v — il or T. 

If we continue the splitting procedure it will give rise to a tree like structure and it is possible to label each sublattice 
by its path in this tree. In other words, we label a sublattice by specifying the sequence of O and T branches that 
constitute the path leading to the sublattice. The idea of above labeling comes from the fact that this branching 
process actually gives a family classification for the sites of the PD chain, and we can consider this tree as the 
genealogical tree for this lattice. So finally we can write GAr(q, cu) as 

G N (q,w) = ^2 P( P ath)G^ ath) {q,uj), (9) 

all paths 

where the sum is over all possible paths in the genealogical tree for a given number of branching. We terminate each 
path in the genealogical tree using the criteria that the corresponding renormalized coupling constants k^ and k$ flow 
to zero at this stage of iteration. In this limit the computation of each term in the Eq. (H) becomes trivial (see Eq. 
(|J)) and one can express every G^""*' 1 into the following general form 

G% ath) - - {x a F*Je* a + x p F* /el + x^/C,) , (10) 

where e*'s and F*'s represent the appropriate renormalized parameters, and x a ,xp and x 7 are the concentrations of 
a, (3 and 7 sites in the PD chain. From Gjv(?,w) the dynamic structure factor can be obtained from Eq. ([!]). 

The merit of this scheme is that one has to generate all possible paths in the genealogical tree using a simple 
algorithm, then iterate the recursion relations To or Tr sequentially along these paths, and finally determine S(q, u>) 
using equations (Q), (^) and There is no approximation involved in this method and the dynamic structure 

factor for quasiperiodic lattices can be obtained with arbitrary accuracy, the accuracy level being set by the smallncss 
of the renormalized values of Ul and k$ ■ 

The dispersion relations for the system can be easily obtained from S(q,u>) using the fact that the lu and q values 
corresponding to every non-zero values of S(q,w) constitutes a point in the dispersion curve. Thus scanning the entire 
lu — q plane for non-zero S(q, to) one can get the dispersion curve for the PD chain, and this gives the whole spectrum 
of normal mode frequencies for the system. 



IV. LATTICE MODES AND THE TRANSFER MATRIX METHOD 



In this section we describe briefly the transfer matrix method for determining the normal mode vibrational states 
in a one-dimensional lattice. The equations of motion for classical vibration of the chain are 

m l w 2 M l = ki,i+i(iH - Ui+i) + ki,i-i(ui ~ Ui-i), (11) 

where Ui is the displacement of the zth atom. In the transfer matrix method one can cast the above equation in the 
following matrix form, 

\ ^ ) \ 1 
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In any normal mode, the displacement u n of any arbitrary atom n can be obtained from the product of a series of 
transfer matrices M„_iM„_2 ■ ■ • MiM\ by specifying the initial displacements u\ and uq- However, in this technique 
the determination of the amplitudes at various sites requires a prior knowledge of the normal mode frequency. Since 
the study of the structure factor gives the entire normal mode frequency spectrum, we can characterize the localized 
or delocalized nature of the lattice modes in quasiperiodic systems from the asymptotic behavior of the amplitudes 
at large distances. 

At this stage, it is worthwhile to mention that in the PD chain a-type sites always occur in pairs, and an analysis 
of this dimmer-type correlation provides an understanding about the extended nature of the electronic eigenfunctions 
in quasiperiodic and disordered systems [p|Jl^|. In a recent work, Dominguez-Adame et al. |l4| have shown that 
this dimmer-type correlation determines the special frequency at which delocalized vibrational mode can exist even 
in a disordered one-dimensional chain. For a PD chain there are three types of transfer matrices M a ,M@ and M 7 
corresponding to the a, (3 and 7 sites in the lattice and it is easy to show that for e a = 0, apart from a constant 
phase factor, the total transfer matrix across the PD chain simply reduces to that of a perfectly periodic chain with 
alternating (3 and 7 atoms. Consequently the system supports delocalized vibrational mode for the frequency lu which 
satisfies the equation e a — 0. Since this a — a pairing is present in all the subsequent renormalized lattices, in principle 
one can determine the full frequency spectrum for delocalized modes by solving the equations e a = e a = e a ■ ■ ■ = 0. 
However, as pointed out earlier, the renormalized e Q 's are polynomials in to, and the degree of the polynomial increases 
very rapidly with the progress of iteration. So it becomes practically impossible to solve these high degree equations 
with adequate accuracy. Thus dimmer approach turns out to be inefficient for obtaining the frequency spectrum of 
the delocalized modes. On the other hand, the dynamic structure factor gives the entire spectrum of normal mode 
frequencies without much effort and we can easily determine the nature of the normal modes using transfer matrix 
method. 



V. RESULTS AND DISCUSSIONS 



With the above background we shall now present in this section the results for some specific cases of the PD chain. 
We begin by considering the on-site model of the PD chain. In the on-site model two types of masses itia and mj, 
connected by identical spring constants, are distributed on a lattice following PD ordering, and the spacing between 
the atoms are equal. We can recover the on-site model by setting m a = m 1 = rriA, mp = rriB, — kg — k and 
o-l = as = a. For numerical calculations, we choose niA = 2,tyib = 1, fc = 1 and a = 1. In the first column of Fig. 
2, we have plotted the dynamic structure factor S(q,u>) as a function of q for three normal mode frequencies, while 
the amplitudes of the atoms in these normal modes are shown in the second column of the figure. The first, second 
and third rows of Fig. 2 respectively correspond to the frequencies io\ — \^lu\ = 0.38196 and lu\ = 1.54739 and these 
frequencies respectively satisfy the conditions e a = 0, e a = and e a = 0. Though these normal mode frequencies can 
be obtained directly from the dynamic structure factor itself, here we have determined them from the above conditions 
for getting better numerical accuracy. This figure clearly depicts that for the on-site model, the dynamic structure 
factor S(q, lu) becomes a periodic function in wavevector q. From Fig. 2a, 2c and 2e we see that the peaks in S(q, lu) 
corresponding to the frequencies ui, lu-i and L03 occur respectively at q = (2n + l)n/2, (2n + l)7r/4 and (2n + l)7r/8, 
where n = 0, 1, 2 • • •. Also the corresponding normal modes display periodicity of periods 4, 8 and 16 in units of the 
lattice spacing. 

Now we shall discuss quite elaborately the underlying physical reason for this periodicity in the dynamic structure 
factor of the PD chain. For this purpose we restrict ourselves only to the frequency lu — lu\. The frequency lu\ is a 
root of the equation e a = 0, and as mentioned earlier, for this frequency one can effectively map the PD chain into an 
ordered binary chain composed of alternating (3 and 7 type atoms. The dynamic structure factor and the dispersion 
relation for this periodic lattice are respectively given by 



S(q,Lu) = — Im 

and 



ep + e 7 + 2(k L + k s ) cos qa 



2(e/3e 7 — k L — k s — 2khks cos 2qa) 



(13) 



lu 1 = (l/2m /3 TO 7 ){(fc L + k s )(m fj + m 7 ) ± [{k L + k s ) 2 {mp - m 7 ) 2 

+ 4m )3 m 7 (fc| / + kg + 2k L k s cos 2qa)} 1/2 }. (14) 

Setting lu = u>x, mp = I, m 1 = 2 and k^/kg = 1 in Eq. (|l3|), in order to afford a comparison with the above 
results presented in Fig. 2, we see that the positions of the peaks in S(q,Lu) for this periodic binary chain are situated 
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at q = [2n + l)ir/2. These positions of the peaks are exactly the same as those observed in Fig. 2a for the on- 
site model of the chain, though their relative intensities may differ. In other words, at w = lo\, the the peaks of 
S(q, co) for the PD chain exactly coincides with that of a periodic lattice whose unit cell contains two atoms (3 and 7. 
From this correspondence we can say that the quasiperiodic PD chain supports delocalized vibrational mode for this 
frequency and the oscillation of atomic displacements in Fig. 2b actually reflects this extended nature of the mode. 
In a similar fashion the periodicity of S(q, co) in Fig. 2c and Fig. 2e can be explained, and we see that for these cases 
the dynamic structure factor resembles that of two other periodic lattices whose unit cells consist of four and eight 
atoms respectively. Thus as far as the dynamics are concerned, quasiperiodic PD chain effectively behaves like an 
assembly of periodic lattices with periods 4, 8, ■ • ■ etc. and it is quite natural that the corresponding normal modes 
are extended. In fact, we have observed that for the present choice of the parameters of the system, every normal 
mode of the PD chain exactly becomes identical to that of an appropriate periodic lattice. 

The above equivalence between the PD chain and a set of periodic lattices become much more apparent from the 
dispersion relation for the PD chain. We have determined the dispersion relation of the PD chain by scanning the 
entire lo ~ q plane for non-zero values of S(q, 00) and it has been displayed in Fig. 3. In this figure we have also plotted 
the dispersion curve (solid lines) of the periodic lattice consisting of two type of atoms (3 and 7. In the periodic case 
with alternating (3 and 7 type atoms, we have one acoustic and one optical branch. Fig. 3 shows that there is an 
acoustic branch in the PD chain which coincides with the acoustic branch of this particular periodic lattice, and both 
of these lattices have a common Brillouin Zone boundary at q = n/2. This indicates that there is an acoustic branch 
of the PD chain which comes from an effective periodic lattice of unit cell size two in terms of the lattice spacing. 
However, in the case of the PD chain, we observe that there are an infinite number of acoustic branches each having 
its own characteristic Brillouin Zone, and Fig. 3 clearly shows that the edges of these Brillouin Zones are situated at 
q = 7r/4, 7r/8, • • ■ etc. In Fig. 3 we have displayed the dispersion curve in the extended Brillouin Zone scheme and 
Zone boundaries other than those mentioned above correspond to the repetition of the various basic Brillouin Zones. 
The overlap of the various acoustic branches implies that these Brillouin Zones have a nested structure. Thus it is 
quite interesting that the quasiperiodic PD chain plays the role of an infinite number of periodic lattices whose unit 
cell sizes increase as 2 n , where n — 1, 2, • • ■ 00. The optical branches of the PD chain are in fact the superposition of 
all the optical branches of these periodic lattices and naturally it would be quite complex. From our calculations we 
really observe that the optical branches have a very complex structure (see Fig. 3) and actually there can be infinite 
number of the optical branches for the PD lattice. We observe the following additional features in the optical branch 
of the PD lattice. For PD lattice we find a well-defined structure in the region of the oj — q plane corresponding to 
the gap of the periodic (3 — 7 binary chain, which is characteristic of the quasiperiodic order of the chain. Fig. 3 also 
shows that such a structure is absent for low values of q. This is due to the fact that in this long wave length limit 
(wavelength A >> a^ or as) the detailed structure of the lattice is not very important, and both the periodic [3 — 7 
chain and the quasiperiodic PD chain exhibit similar behavior. 

In Fig. 4 we have displayed S(q 1 uo) as a function of q for other models of the PD chain including some additional 
graphs corresponding to the on-site model. The graphs for the on-site model are shown in Fig. 4a and 4b corresponding 
to the frequencies uj 2 = 2.26829 and oo 2 — 1.65555 respectively. The graphs for the bond-model are shown in Fig. 
4c and 4d which has been realized by setting m a = mp = m 7 , ^ ks and 7^ as- In this case we choose 
m a = mp = m 7 = 1, ki, = 1; ks — 2, — 2 and as = 1 and the Fig. 4c and 4d respectively correspond to the 
frequencies oj 2 — 2.0 and lo 2 = 3.09447. In Fig. 4e and 4f we have also presented the results for the most general case 
of the PD chain and here the parameters are taken as m a = 1, mp = 2, m 7 = 3 (all in arbitrary units), ks/kL = 2 
and a^/as = 2. The graphs in Fig. 4e and 4f are plotted for frequencies lo 2 = 2.0 and u 2 = 2.42427 respectively. The 
most striking feature is that in every case we observe that S(q, lo) is always periodic in wave vector q. So we can say 
that this periodicity is a characteristic feature of the PD lattice and it is quite independent of the specific choice of 
the parameters of the system. 

We have demonstrated the underlying physical reason for this periodicity in S(q, to) only for the on-site model of 
the PD chain, which owes its origin to the fact that for every normal mode frequency one can exactly map the system 
to various periodic lattices. Similar kind of analysis is also possible in other cases, the unit cell structure being quite 
complex in these situations. The most interesting consequence of the periodicity of S(q,co) with respect to q is that 
the system supports delocalized vibrational modes even in absence of the translational invariance. 

At this stage we would like to make the following remark. It should be noted that there exists a polynomial invariant 
/ associated with the recursion relations Eq. (^) and (||) for the PD chain which is given by [|| 

/ = {(e a - e 7 )(e Q - ep) -k%- k 2 )/{2k L k s ) + 1. (15) 

Depending upon the choice of parameters of the system, / may be independent of lo as well as it may be a function of 
to. In the electronic case, one normally equates the invariant to zero for finding the energy eigenvalues of the extended 
states of quasiperiodic systems, provided those energies are allowed ones. Similarly we can obtain the delocalized 
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normal mode frequencies from the condition I — 0. However, we have to check whether the frequencies obtained from 
the above condition give allowed eigenmodes or not. For m a = l,mp = 2,m 7 = 3, ks/k^ = 2 and a^/as = 2 it 
is found that the condition I = is satisfied for lo 2 = 1.5. But this is not an allowed frequency of the system (see 
Ref. which is also confirmed from the calculation of S(q,u) and it turns out that it is always zero for every q 

in this case. However, depending on the choice of the parameters of the system, the condition 1 = may give a few 
allowed delocalized normal mode frequencies. This condition actually corresponds to the commutation of the transfer 
matrices for two successive generations of the chain, and thus leads to delocalized states. It should be noted that the 
commutation conditions for various successive generation always give the same set of frequencies due to the presence 
of the invariant in this system (see Ref. [Q for details). Our RG analysis, however, shows that an infinite number of 
delocalized eigenmodes are possible in the PD chain due to a — a pairing at all length scales. 



VI. CONCLUSIONS 



In conclusion, we should stress that the renormalization group method is the right tool for taking into account the 
symmetry of quasiperiodic systems and the present RSRG scheme offers a very efficient method for computing the 
dynamic structure factor of these systems. The merit of the present scheme is that it provides an exact method for the 
determination of S(q, u>) of infinite self-similar quasiperiodic lattices, and in practice, one can calculate S(q, u>) with 
arbitrary accuracy. This method essentially involves the iteration of certain recursion relations, and so it enormously 
reduces the computational task. We have applied this method to the quasiperiodic PD chain and obtained several 
new results. Apart from computing the dynamic structure factor, we have also determined the dispersion relation of 
this lattice. The most interesting result is the observation of periodicity in S(q, u>) with respect to q for every normal 
mode frequency, and this periodicity actually gives the signature for the presence of normal mode vibrations which 
are delocalized over the entire lattice. A direct consequence of the periodicity of S(q, u>) is that one can infer about 
the existence of the delocalized vibrational modes in the system simply by looking for the periodicity in the scattering 
intensity data from diffraction experiment measurements. Another important outcome of the present work is that the 
entire spectrum of frequencies for delocalized modes can be obtained very easily from the dynamic structure factor, 
whereas the determination of the full spectrum by the method of Ref. |Q is practically impossible. 
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FIG. 1. A section of the period-doubling chain illustrating the sublattice splitting. The symbols •, A and O respectively 
represent a, (3 and 7 sites. 

FIG. 2. Plot of S(q, lj) versus q, and u„ versus n for the on-site model with m a = m 7 = 2, mp = 1 and k^jks = 1, a^/as = 1. 
First, second and third rows respectively correspond to u> 2 — 1.0, lu 2 — 0.38196 and lo 2 — 1.54739. 

FIG. 3. The dispersion relation for the period-doubling lattice with the same parameters as those in Fig. 2. 

FIG. 4. Plot of S(q,u>) versus q. On-site model with parameters ra a — m 7 = 2, mp — 1, kh/ks = 1 and cll/cls = 1 for (a) 
lo 2 — 2.26829 and (b) lj 2 — 1.65555. Bond-model with parameters m a — rap — m 7 = 1, ks/kh = 2 and ah/as = 2 for (c) 
uj 2 = 2 and (d) lu 2 — 3.09447. Mixed model with parameters m a — 1, rap — 2, m 7 = 3, k$/k L — 2 and a L /as = 2 for (e) 
lj 2 = 2 and (f) uj 2 = 2.42427. 
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